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Results of recent density functional theory calculations 
for alkali metal adsorbates on close-packed metal surfaces 
are discussed. Single adatoms on the (111) surface of Al 
and Cu are studied with the self-consistent surface Green- 
function method by which the pure adsorbate-substrate in- 
teraction may be analyzed. Higher coverage ordered adlay- 
ers of K on Al(lll), Na on Al(lll), and Na on Al(OOl) 
are treated using the ab-imtio pseudopotential plane wave 
method which affords the prediction of coverage dependent 
stable and metastable adsorbate geometries and phase tran- 
sitions of the adsorbate layers. Together, these studies give in- 
sight and understanding into current key issues in alkali metal 
adsorption, namely, the nature of the adsorbate-substrate 
bond at low coverage and the occurrence of hitherto unan- 
ticipated adsorbate geometries, and the associated electronic 
properties. 



I. INTRODUCTION 

Alkali metal adsorption systems have served as a theo- 
retical paradigm for understanding the properties of ad- 
sorption on crystaUine sur faces 0-0. This is due to the 
simple electronic structure of the alkali metal atoms and 
also because it was understood that the adatom would 
not complicate the system of interest by intermixing with 
the substrate, but would stay on the surface and occupy 
a high-symmetry, highly coordinated position. Early ex- 
periments revealed that these "simple" adsorbates had a 
marked effect on the properties of the surface on which 
they were adsorbed [0. For example, in the seminal 
work of Taylor and Langmuir jlj] for the Cs/W system, 
it was shown that the work function of the clean sur- 
face is reduced by several electron volts on Cs adsorption 
and that at low coverage almost all the Cs adatoms are 
desorbed as positive ions. In the light of such results, 
the alkali metal bond to the substrate was viewed as a 
spectacular example of ionic bonding akin to that in al- 
kali halides . As a consequence of these distinct and 
general properties, alkali metal adsorption has come to 
play an important role in technological applications. For 
example, the desirable properties of the considerable en- 
hancement of the electron emission of cathodes and the 
improved activity and selectivity of heterogeneous cata- 
lysts (ijjl^ are both achieved by the addition of alkali 



metal atoms to surfaces |16 



A simple picture of the interaction between alkali 
metal atoms and a metal surface, and of the resulting 
chemisorption bond was proposed by Langmuir ||l7|| . He 
assumed that the alkali metal atom transfers completely 
its valence s electron to the substrate. In a more rigorous 



theoretical description of alkali metal adsorption Gur- 
ney Q, in 1935, proposed a quantum mechanical picture 
applicable at low coverages, where the discrete s level 
of the free alkali metal atom broadens (and shifts) be- 
coming partially emptied as a result of the interaction 
with the substrate states as it approaches the surface 
(see Fig. As follows from this picture the alkali metal 
atom is partially positively charged which induces a neg- 
ative charge density in the substrate giving rise to an ad- 
sorbate induced dipole moment which naturally explains 
the reduction in the work function of clean surfaces. The 
building up of such adsorbate induced dipole moments 
also leads to the understanding that the dominant in- 
teraction between the adsorbates is repulsive and that 
the adatoms should form a structure with the largest 
possible interatomic distances compatible with the cov- 
erage. In this description it is expected that with in- 
creasing coverage the adsorbate-adsorbate distance grad- 
ually decreases, and the electrostatic repulsion between 
the adatoms increases. To weaken this repulsion, i.e., 
to lower the total energy, some fraction of the valence 
electrons flow back from the Fermi level of the metal 
to the adsorbate. Thus, a reduction of the adsorbate- 
induccd dipole moment, i.e., depolarization takes place. 
This picture of Langmuir and Gurney has been held with- 
out question for many years. 

The Langmuir-Gurney model was fully supported by 
systematic studies of Lang and Williams , who per- 
formed density functional theory calculations of isolated 
adatoms bonded to the surface of a simple metal ("jel- 
lium"), and by Lang's earlier work on the alkali 
metal induced work function change. In these calcula- 
tions the ions of the metal substrate were approximated 
by a semi-infinite uniform positive background which ter- 
minated abruptly along a plane, otherwise known as the 
"jellium" model. A jellium model was expected to rep- 
resent well the electron density of a free-electron-like s- 
p bonded metal [Q, and indeed these studies demon- 
strated that many ground-state properties of metal sur- 
faces and metal-adatom systems can be described in a 
physically transparent manner. As an example, we show 
in Fig. Ijthe change in work function with increasing cov- 
erage for parameters corresponding to Na on Al(lll). 
The curve possesses a form similar to that often ob- 
served experimentally. This "usual" shape of A$(8) is 
explained as a consequence of the above-mentioned depo- 
larization of the alkali metal induced surface dipole mo- 
ment induced by continuous reduction of the adsorbate- 
adsorbate distance and corresponds to a rapidly decreas- 
ing work function at low coverage, reaching a minimum 
at about Q = 0.15, and subsequently rising towards the 
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value of the pure alkali metal. 

The above description represents a sketch of the "tra- 
ditional view" of alkali-metal adsorption, which, for ex- 
ample, is underlying the analysis of experimental re- 
sults in various articles in the book of Bonzel, Bradshaw 
and Ertl However, since 1991 it has become clear 

through density functional theory studies |^,|lj , as well 
as through several experimental investigations, that the 
interactions and reactions of alkali metal atoms at metal 
surfaces are more complicated and exciting than hitherto 
expected. The calculations, which take the full atomic 
structure into account including relaxation of atoms from 
their ideal positions as well as surface reconstructions, 
show that the traditional view of alkali metal adsorption 
outlined above, is only part of the whole picture, and that 
we must now adopt an updated view. Phenomena such 
as the following may occur: 1) The alkali metal atoms 
may not necessarily assume highly coordinated sites on 
the surface ||2l|-|3^. 2) The alkali metal adatom may kick 
out surface substrate atoms and adsorb substitutionally 
pO|,pl 30-3^. Substitutional adsorption has been shown 
in three cases to occur as the result of an irreversible 
phase transition from an "on-surface" site by warming 
to room temperature, without change in the periodicity 
of the surface unit cell or of the coverage |3^-|3^] . 3) The 
alkali metal atom may switch site on variation of cov- 
erage [E9l^MQ-M and 4) island formation may occur 
Ill3| , t3l|j44[ - ^ ]. 5) There may be a strong intermixing of 
the alkali metal atom with the substrate surface [p9|-p2[ , 
as for example the formation of a four layer ordered sur- 
face alloy which had been identified for Na on Al (111) at 
a coverage 0Na = 0.5 |^,|2|. The coverage is defined in 
this paper such that for = 1, the number of adatoms 
is the same as the number of atoms in the clean unre- 
constructed surface layer. Alkali metal induced surface 
reconstructions are well known on the more open surfaces 
1^-^ as these clean surfaces are close to a structural 
instability [^6|. But that significant reconstruction can 
occur on close-packed fee (111) and (001) surfaces had 
not been expected previously. 

The associated electronic properties of these hitherto 
unexpected surface atomic arrangements, perhaps not 
surprisingly, also exhibit behavior deviating from expec- 
tations based on early ideas. For example, experimental 
measurements of the change in work function as a func- 
tion of alkali metal coverage can be quite different to the 
"usual" form of Fig. ||, for example, they may not ex- 
hibit a minimum at all, or may have a minimum at a 
significantly higher coverage. In the same manner, the 
density of states induced by alkali metal adsorbates may 
not always correspond to that expected from the model 
of Gurney. 

Contrary to the idea mentioned above that the tra- 
ditional picture be kept as representing part of a re- 
vised view, are recent theoretical studies which sug- 
gested a different picture of alkali metal adsorption be 
adopted | |57| , p8[ . In this picture there is practically no 
charge transfer from the alkali metal adatom to the sub- 



strate but rather a polarization of the adsorbate and the 
adsorbate-substrate bond is classified as covalent as op- 
posed to the general and widespread understanding that 
it is of ionic nature at low coverages. This interpretation 
has seemingly been supported by experimental studies 
to the extent that it was concluded that the traditional 
picture be entirely discarded p9[ . This is a controversial 
issue but we consider it a largely artificial conflict based 
on semantic differences. For a recent discussion we refer 
to Ref. | ]60| , |6l| , and a short summary of some of these 
arguments is given in Section III below. 

The purpose of the present paper is to discuss results 
of recent ab-initio calculations for alkali metal adsorbates 
on close-packed metal surfaces. These density functional 
theory studies provide an accurate description of the mi- 
croscopic electronic structure of the adatom as well as of 
the atomic geometry. In Section II we give a brief out- 
line of the theoretical methods used. The results of Lang 
and Williams [^,D, Bormet et al. j6l), and Yang et al. 
p2[ , for calculations of single adsorbates on jellium, on 
Al(lll), and on Cu(lll), respectively are discussed in 
Section III. These calculations afford an analysis of the 
nature of the adsorbate-substrate bond without it being 
obscured by the influence of other adsorbates. The re- 
sults of these studies do not support a dismissal of the 
traditional view and the adoption of a covalent picture of 
alkali metal adsorption when adsorbed as isolated atoms. 
Rather, the results show that the traditional view is valid 
and useful but should only be applied in the limit of 
very low coverage and non-activated adsorption (see Sec- 
tion IV for the definition of activated and non-activated 
adsorption). The comparison of the results for different 
substrates enables one to establish the influence of the 
atomic structure and chemical nature of the substrate on 
the calculated quantities. In Section IV we discuss recent 
results obtained using the self-consistent pseudopotential 
plane wave method for higher coverages of alkali metal 
adsorbates on Al. From these calculations "unusual" ad- 
sorption sites and phase transitions have been predicted. 
The calculations also provide insight and understanding 
into the reasons why the unexpected adsorbate arrange- 
ments occur. In subsections IV A and IV B we discuss 
respectively, the coverage dependence of the adsorbate 
geometry, and the adsorbate-adsorbate and adsorbate- 
substrate bonding and the associated electronic proper- 
ties. Section V contains the conclusion. 



II. THEORETICAL METHODS 

To provide some understanding of the two different 
calculation procedures, namely, the self-consistent sur- 
face Green-function (SSGF) method and the pseudopo- 
tential plane wave method, used to obtain the results dis- 
cussed in the present paper, we give a short description 
of the theoretical basis. Both methods employ density 
functional theory (DFT) with the local density approx- 
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imation (LDA) for the exchange-correlation functional 
(DFT-LDA). We will first give an outline of this, then 
explain each of the calculation procedures separately. In 
the equations below Hartree atomic units are used. 



A. Density functional theory 

Density functional theory, as introduced by Hohen- 
berg, Kohn, and Sham []63|-|65[ has proved to be a con- 
ceptually and practically useful method for studying the 
electronic and atomic structure, as well as elastic and 
vibrational properties of many-electron poly-atomic sys- 
tems. This theory represents a firm and exact theoretical 
foundation whereby all aspects of the electronic struc- 
ture of a system in a ground state are completely deter- 
mined by its electron density n{r). The total energy of 
the ground state is given by 

E^°^[n{r)] = T[n{r)] + [ V{r)n{r)dr 



and 



n(r)n(r') 



(1) 



where V{r) is the "external potential" due to the atomic 
nuclei (or the ion cores) and T[n] is the kinetic energy 
functional for non-interacting electrons having the den- 
sity n(r). The functional i?a;c[?i] contains all the quantum 
mechanical many-body effects and is called the exchange- 
correlation energy. Both T[n\ and i?2;c[?^] are not known 
explicitly but at least the kinetic energy can be evalu- 
ated exactly when the density is written as a sum over 
orthonormal single-particle functions, 



N 



n(r) = 



(2) 



1=1 



where N is the number of electrons in the system. This 
gives, 



N 



(3) 



The variational property of the total energy then leads 
to the equation. 
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where 



V,a{r)^V{r) + VHir) + V^c{r) 



(4) 



(5) 



which is to be solved self-consistently. The Hartree and 
exchange correlation potentials are. 



dn{r) 



The Coulomb energy associated with interactions among 
the M nuclei (or ions) at positions R/ is given by 



M,AI 

E 



ZrZj 



(6) 



The only unknown quantity is the exchange-correlation 
energy functional, Exc[n\, and the quality of the solution 
of the full many-body problem is only limited by the 
quality of the approximation used for it. The exchange- 
correlation energy can be written as, 



Exc[n{r)] = J exc[n{r)] n{r) dr 



(7) 



A simple yet effective way to evaluate Exc[n] is to use 
the LDA which means that the exchange-correlation 
energy density functional, e^c [n] , is approximated by the 
corresponding expression for the homogeneous electron 
gas. That is, 



,LDA 



in) 



(8) 



where e^Jp^('T-) is the exchange-correlation energy per 
particle of the many-body system which is known most 
accurately, namely, the uniform electron gas. Expe- 
rience from numerous calculations show that for well 
bonded situations the LDA represents a good descrip- 
tion of the quantum-mechanical many-body interactions 
in poly-atomic systems, and calculated electron densities, 
n(r), atomic geometries, and total energy differences are 
in fact very reliable. In actual calculations, often some 
additional, more severe numerical approximations are ap- 
plied. The latter can (and should) be carefully tested. 



B. The self-consistent surface Green-function 
method 

The self-consistent Green-function (SSGF) method 
is outlined briefly here and is presented in detail in 
Refs. [ ]60| , plp^ ]. It gives an accurate and efficient descrip- 
tion for isolated adsorbates on semi-infinite substrates. 
Using this method the total and adsorbate induced den- 
sity of states (DOS), charge densities and density dif- 
ferences can be obtained, as well as the total energy, 
total energy change, and forces on the adsorbate. In 
this approach the system is split into a reference sys- 
tem and a perturbative part with a perturbation poten- 
tial, A T4ffM, which has to be evaluated self-consistently. 
The reference system is taken to be the ideal, clean 
surface and the corresponding Green- function, G", is 
calculated employing the layer-Korringa-Kohn-Rostoker 
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(KKR) Green-function method |67|. The Ml Green- 
function of the adsorbate system is given by the Dyson 
equation, 



G G° + G° AVcff G 



(9) 



Subtracting G° from both sides of equ. ^ gives, 

AG = G - G° = G° AV^ff G (10) 



= G°Ay(l-G°Ayeff)"^G° 



(11) 



The perturbation potential, Al^g , is written as, 

A^effKr)] =Foff[n(r)]-yO , (12) 

where V'^ is the potential defining the reference system 
and the Green function G°, and n(r) is the density of 
the valence electrons. Thus, V^ff[n(r)] is the functional 
of the effective potential as described by the Kohn-Sham 
equation (equs. ^ and |^). The perturbation potential, 
A14fi, is strongly localized in real space. This localiza- 
tion allows equ. |ll] to be solved efhciently in a rather 
small, localized basis set. In the work discussed below, 
s,p, and d gaussians and three different decay constants 
per angular momentum number are used which are cen- 
tered on the adatom and on the three nearest-neighbor 
substrate atoms 61, 6|]. 



Once the change in the valence electron density is cal- 
culated self-consistently, the total energy, , as well as 
the forces, F 



on the adsorbate as a function of 



its position R can be evaluated |61| 



C. The ab-imtio pseudopotential plane wave method 

A very successful way to solve the Kohn-Sham sin- 
gle particle equation and to perform accurate calcula- 
tions for a periodic system is achieved when ion cores 
are treated by ab-initio norm-conserving pseudopoten- 
tials and a plane waves basis set is employed. The great 
advantage of this approach is that the basis set is unbi- 
ased, i.e., the numerical accuracy can be tested system- 
atically. Furthermore, it provides accurate forces on the 
atoms without complications typically encountered when 
other basis functions are used. The wave functions are 
written as. 



■!/'j(k,r) = ^ Cj,G+kexp(i(G -t- k) • r) 



(13) 



G+k|<G„ 



where G are the reciprocal lattice vectors. Because the 
coefficients, Cj^c-i-k, for plane waves with high kinetic 
energy decrease and become negligible, the plane wave 
basis set can be truncated to include only plane waves 
that have kinetic energies less than some particular cut- 
off energy, ^|GmaxP = Ecut- Thus, the cutoff energy 
defines the quality of the basis set and it is one of the 



important numerical parameters of the calculations for 
which it must be carefully tested that they do not affect 
the results. In the pseudopotential plane wave method 
a crystal surface is modeled by a large supercell which 
contains a slab a certain number of atomic layers thick 
which is surrounded by vacuum. This supercell is re- 
peated throughout space. The surface is then periodic in 
the plane of the surface on the scale of the supercell, but 
the periodicity perpendicular to the surface, imposed by 
the supercell, has no influence because for a sufficiently 
thick vacuum region the different slabs will not inter- 
act. However, care must be taken that the slab is thick 
enough so that its two surfaces are practically decoupled 
from each other. For a further description of the numer- 
ical details we refer to Refs. mJGSI. 



III. ISOLATED ADSORBATES 

By studying isolated adatoms on surfaces the nature of 
the pure adsorbatc-substrate bond can be analyzed. Self- 
consistent DFT-LDA calculations for isolated adatoms 
on a semi- infinite jellium surface have been performed 
by Lang and Williams ||,^. Bormet et al. carried 
out an analogous study for a representative set of iso- 
lated adsorbates on Al (111). In the latter work, instead 
of using a jellium model for the Al substrate, the atomic 
structure of the substrate is taken fully into account. 
Therefore, any significant differences in the results of the 
two studies is directly attributable to the treatment of 
the substrate. As a complementary study to the adsor- 
bate on Al(lll) systems, Yang et al. investigated 
isolated adsorbates on Cu(lll). From comparison with 
the results for Al(lll), this work allows the influence 
of the Cu id electrons on the adsorbate bonding to be 
identified. In both of these latter calculations the self- 
consistent surface Green-function (SSGF) method was 
used. The set of adsorbates from the third row of the 
periodic table, were chosen in order to establish chemical 
trends. We concentrate here, in particular, on the case of 
alkali metal adsorption, but briefly describe the results 
for the other adsorbates, namely, Si and CI in order to 
highlight the very different and characteristic behavior 
of the adsorbates which is found to be highly correlated 
with the electronegativity of the adsorbate with respect 
to the substrate. 



A. Nature of the bond in the limit of very low 
coverage 

For adsorption on Al (111) and Cu (111) the single ad- 
sorbates were placed in the fcc-hollow site, and in or- 
der to compare with the results of Lang and Williams 
|^,||, the surface was taken to be unrelaxed. This is, in 
fact, a good approximation because for non-activated ad- 
sorption the actual substrate relaxation is rather small. 
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The vertical height of the adsorbates was optimized by 
searching for the minimum of the total energy. From 
the study of Ref. it was found that the equilibrium 
heights of the adsorbates above the surface are smaller 
for the Al (111) substrate than for the jcUium substrate. 
With respect to the center of the top aluminum layer 
the heights are approximately 3, 20, and 20% smaller on 
Al (111) than on jellium for Na, Si, and CI, respectively. 
For jellium the "center of the top layer" is defined to 
lie half a bulk inter layer spacing of Al(lll) below the 
jellium edge. The larger values obtained in the jellium 
calculations are due to neglecting the atomic structure of 
the substrate which has the ability to bind the adsorbates 
more strongly. It is interesting that for Si on a jellium 
surface and using first-order perturbation theory for the 
crystal lattice, an equilibrium height is obtained which is 
identical to the full SSGF result H,^. On comparison 
of the results for adsorption on the Cu and Al substrates 
it is found that the jellium— ^aluminum trend continues: 
the effective radius of the studied adsorbates on Cu (111) 
is smaller than on Al (111) by approximately 11, 15, and 
17%, respectively, for Na, Si, and CI. Clearly, in this case 
it is more meaningful to compare effective radii than in- 
ter layer spacings due to the different lattice constants of 
Al and Cu. The smaller values reflect the existence of a 
stronger bonding of the adsorbates to the substrate for 
Cu(lll) as compared to Al (111). This is identified as 
being due to the Cu d electrons |62| . 

Adsorbate induced density of states (DOS) and charge 
density difference plots are informative in discerning the 
nature of the adsorbate-substrate bond. The adsorbate 
induced DOS indicates whether the state of interest is 
occupied, unoccupied, or partially occupied, while charge 
density difference plots show whether there is a shift of 
charge density from the adsorbate towards the substrate, 
vise versa, or a sharing of charge between the adsorbate 
and the substrate. 

The adsorbate induced DOS obtained from the SSGF 
calculations of Bormet et al. j6l) and those of Lang and 
Williams |^ for adsorbates on jellium, corresponding to 
a high electron density metallic substrate (rg — 2 bohr), 
are shown in Figs. ^ and respectively. The adsor- 
bates investigated in the work of Bormet et al. |6l[| were 
Na, Si, and CI and in that of Lang and Williams |9| they 
were Li, Si, and CI. The results of both of these cal- 
culations agree qualitatively and may be interpreted as 
follows: For both Na and Li the adsorbate resonance lies 
well above the Fermi level and is thus largely unoccupied. 
This indicates that the valence electron of the alkali metal 
atom (or part thereof) has been transferred to the sub- 
strate and the adatom is partially positively charged. In 
an opposite manner, on adsorption of chlorine, the reso- 
nance in the curve corresponding to the CI 3p resonance 
lies well below the Fermi level. Thus, this resonance is 
fully occupied which implies that a transfer of electronic 
charge density from the substrate to the CI adatom has 
taken place; the adsorbed CI atom is partially negatively 
charged. For the adsorption of an isolated Si atom it can 



be seen from the jellium calculations that the Si 3p reso- 
nance lies just at the Fermi level, which implies that it is 
about half occupied. Lang and Williams |^ have shown 
that the states of the energetically lower half of the p- 
resonance are bonding between the adatom and the sub- 
strate and that the energetically higher states are anti- 
bonding. Because the Fermi level cuts the p-resonance 
approximately at its maximum, the bonding nature of Si 
is covalent. The results for Si on Al (111) show also that 
the bond is covalent. In this case, however, more struc- 
ture occurs in the DOS than in the jellium calculations. 
This arises because the atomic structure of the substrate 
leads to band structure effects (clearly reflected by the 
structure of the bulk DOS at Ep in Fig. 3a), and these 
induce a splitting of the bonding and anti-bonding states, 
and the adatom density of states exhibits a minimum at 
the Fermi level |^l[. Similarly to the jellium substrate, 
the Fermi level cuts the Si 3p induced DOS roughly in 
the middle. 

The behavior identified in the adsorbate induced DOS 
triggers the distribution of the valence electrons. In Fig. ^ 
we show the valence electron density for the chemisorp- 
tion of Na, Si, and CI on Al(lll). The charge transfer 
from the Na atom to the substrate, as reflected by the ad- 
sorbate induced DOS, is clearly visible in Fig. ^ where it 
can be seen that there is an increase of the electron den- 
sity between the Na adatoms and the topmost layer of 
the Al substrate, and that on the vacuum side of the Na 
atom the electron density has almost disappeared. For 
Si, as expected from the adsorbate induced DOS, a di- 
rectional covalent bond is present between the Si adatom 
and the nearest-neighbor substrate atom. Furthermore, 
it can be seen that the maximum of the charge density of 
the chemisorption bond is closer to the more electroneg- 
ative Si atom. In the case of CI on Al the charge density 
distribution around the adsorbate is almost spherical. 

In Fig. H we consider density difference plots. These 
plots display the electron density of the adsorbate sys- 
tem minus the density of the clean substrate and minus 
the electron density of the free adatom. The results for 
the Al and Cu substrates show a number of similar fea- 
tures, and some interesting differences. Firstly, it can be 
noted that in each case the perturbation to the system 
caused by the adsorbates does not reach very far into 
the metal substrate. The interior is essentially identical 
to that of the clean surface for layers deeper than the 
second. We emphasize that this localization holds for 
the electron density perturbation but not for individual 
wave functions. For Na adsorption, perhaps the most 
noticeable feature is the charge cloud between the ad- 
sorbate and the substrate which reaffirms that charge is 
displaced from the vacuum side of the adsorbate towards 
the substrate side. These results are consistent with a 
partially positively charged adatom which sits on a neg- 
atively charged site. This constitutes a surface dipole 
which locally decreases the work function. The opposite 
situation is found for CI, which is negatively charged and 
sits on an adsorption site that is positively charged. In 
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this case charge has moved from the substrate towards 
the CI atoms and the local work function therefore will 
be increased relative to the value of the clean surface. 
For both substrates Si appears covalently bound with a 
slight electron transfer towards the adatom. Comparing 
in more detail the results for Al and Cu, we see that there 
is more charge between the Na adatom and the top layer 
of the Cu substrate. Also, there is a greater depletion 
of charge at the vacuum side. These effects are consis- 
tent with the fact that the electronegativity difference 
between Na and Cu (0.93 — 1.90 = —0.97) is larger than 
that between Na and Al (0.93 - 1.61 = -0.68). The Si 
atom, which clearly forms a covalent bond with both sub- 
strates, is slightly more electronegative than Al (by 0.27). 
In this respect, it can be seen that more charge resides 
on the Si atom when on Al than when on Cu for which 
the electronegativity difference is zero. The result for CI 
on Cu displays some features similar Si on Cu, namely, 
the pile up of charge density between the CI adatom and 
the nearest Cu atom. This structure is practically absent 
for CI on Al. In all cases for the adsorbates studied, the 
adsorption on Cu exhibits some structure in the valence 
electron density change near the nucleus of the Cu atom 
closest to the adsorbate. This reflects the participation 
of the Cu d electrons in the bonding. 

To summarize. Green-function calculations for single 
adatoms on jellium, Al(lll), and Cu(lll) substrates 
show qualitatively the same trend for the different ad- 
sorbates and this is well described by the adatom elec- 
tronegativities compared to the electronegativity of the 
substrate atoms. For single allcali metal atoms on these 
surfaces, all the calculations support the charge transfer 
picture emphasized in the model of Gurney as providing 
a useful and appropriate description of the situation for 
alkali metal adsorption at very low coverage. Admittedly 
(see also Refs. |^,^) the picture is useful and appropri- 
ate, but it is not unique. A (partially) charged entity in 
front of a metal surface will necessarily induce a screen- 
ing charge density at the metal surface - the quantum 
mechanical realization of the image charge. This screen- 
ing charge density is located between the adatom and the 
substrate surface. Therefore, from a purely mathemat- 
ical point of view, the adatom induced electron density 
can be described in three different ways: i) in terms of 
orbitals of the adatom, ii) in terms of orbitals of the sub- 
strate, and in) by taking both orbital sets into account. 
This ambiguity is due to the fact that the sum of the 
two Hilbert spaces of the adatom and of the substrate is 
over-complete. Whereas a mathematical analysis of the 
adsorbate induced electron density leaves this freedom, 
the chemical interpretation of the DOS, as well as its 
effect on the charge density, as discussed above, clearly 
support the qualitative picture of Gurney. We note, how- 
ever, that this picture applies only to low coverages and 
an unreconstructed substrate. Otherwise, as will be seen 
in Section IV, some significant modifications take place. 



IV. ALKALI METAL ADLAYERS ON AL 

One of the interesting characteristics of alkali metal 
adsorption is the drastic change of several physical quan- 
tities as a function of coverage. The most obvious of these 
properties are perhaps the surface atomic geometry and 
the work function. With respect to the former, different 
commensurate phases are formed by alkali metal adsor- 
bates on metal surfaces. Some of these structures can be 
understood as being due to a repulsive lateral adsorbate- 
adsorbate interaction, causing the adatoms to sequen- 
tially arrange into commensurate phases of ever growing 
compactness |l3| . In addition, as mentioned in the intro- 
duction, behavior quite different to these arrangements 
has been reported, as well as a new class of unexpected 
geometries becoming evident, namely, those of substitu- 
tional and surface alloy structures. 

To learn about the strong coverage dependence of 
the surface atomic arrangements and related properties, 
calculations for higher coverages of alkali metal adsor- 
bates on metal surfaces have been performed for K on 
Al (111) [^ MP, Na on Al (111) |l|,|l],|l]j70| , and Na 
on Al(OOl) I^^Q. The calculations show that on in- 
creasing the coverage from very low coverage, where the 
alkali metal adatoms are uniformly spread over the sur- 
face, to higher coverages, the adsorbate-adsorbate inter- 
actions weaken the adsorbate-substrate bond because ad- 
ditional bonds are built up within the adlayer. We note 
that this is a general feature of bonding in poly-atomic 
systems: When the valence electrons are shared among 
many bonds, each single bond becomes weaker, and an 
increase of the bond length results. In the case of alkali 
metal adlayers this implies that a formation of lateral 
bonds in the adlayer will weaken the adsorbate-substrate 
bonds and increase the adsorbate-substrate bond length. 
This interplay/competition between adsorbate-adsorbate 
and adsorbate-substrate bond formation can give rise to 
interesting adsorbate arrangements and structural phase 
transitions. The results for the three systems are com- 
pared which affords changes due to different alkali metal 
atoms on the same surface to be identified as well as 
changes due to different substrate surfaces for the same 
adsorbate. 

It is necessary to distinguish "non-activated" and "ac- 
tivated" adsorption. The former corresponds to adsorp- 
tion on an only relaxed surface (small rearrangements 
of substrate atoms) where no energy barriers need to be 
overcome to reach the adsorption geometry. The latter 
refers to adsorption involving activated processes where 
it is likely that an energy barrier must be overcome in 
order to reach the adsorption geometry, and the sub- 
strate surface undergoes a significant reconstruction. In 
the present work, non-activated and activated adsorption 
may be regarded as corresponding to structures formed 
at low temperature and at room (or higher) temperature, 
respectively. 



6 



A. Adsorbate geometry 

By calculating the adsorption energy for various cov- 
erages and adsorption sites, detailed insight about sta- 
ble and metastable adsorbate geometries and structural 
phase transitions can be obtained. Furthermore, an ac- 
curate knowledge of the adsorbate and substrate atomic 
geometry is crucial for any additional analysis of the ad- 
sorbate properties, as for example, the surface electron 
density of states, chemical reactivity, and adsorbate in- 
duced work function changes. The technical details of the 
den sity functional theory calculations are given elsewhere 

The stable and metastable geometries are determined 
by an automated movement of the atoms of the adsorbate 
and substrate surface on the Born-Oppenheimer energy 
surface into the total energy minimum. For the analy- 
sis, it is necessary to define the adsorption energy per 
adatom. We will write this definition down using in- 
dices which refer to the adsorption of Na on Al(lll), 
but a transformation to other systems is obvious. The 
adsorption energy is the difference of the total energy of 
the adsorbate system and the total energy of the clean, 
unreconstructed Al(lll) substrate together with a free, 
neutral Na atom. Because total energies of atoms have 
negative values, the present and usually applied defini- 
tion, associates exothermal adsorption with a positive 
adsorption energy j?^. The adsorption energy for the 
on-surface (non-activated) adsorption is given by 

^Na/Al(lll) _ ^Al(lll) _|_ ^Na-atom _^Na/Al(lll) 

(14) 

where £;Na/Ai(iii) ^-j^g total energy per adatom of the 
adsorbate system, E^^ '^^^^^ is the total energy of the clean 
Al(lll) sulastrate, and i5;Na-atom jg ^^^g^^ energy of a 
free Na atom. 

For substitutional adsorption the definition is in prin- 
ciple identical, however, one has to take into account that 
the total energy £;Na/Ai(iii) corresponds to a sys- 
tem with a substitutional Na atom for which the kicked 
off surface Al atom has diffused to a kink site at a step 
where it is rebound. Thus, for the substitutional adsorp- 
tion the adsorption energy is 

^Na/Al (lll)-sub _ ^Al(lll) _|_ ^Na-atom _ ^Na/Al (lll)-sub 
- (^^'""*°'" - ^coh) . (15) 

Here we have written the total energy of the displaced Al 
atom which is rebound at a kink site, as the sum of the 
total energy of the free Al atom and the negative of the 
cohesive energy . Thus, we assume thermal equilibrium 
and take into account that this is transmitted by kink 
sites at steps. The quantity £;Na/Ai(iii)-sub || -g 

the total energy of the slab with the adatoms adsorbed 
in substitutional sites. 



It is sometimes useful for a more transparent discus- 
sion to describe the substitutional adsorption in terms of 
two independent processes, namely, at first the creation 
of the surface vacancy structure, and then the chemisorp- 
tion of the adatoms on the substrate surface with vacan- 
cies. Firstly, the energy required to create the surface va- 
cancy structure is equal to the total energy of the (recon- 
structed) substrate (with surface vacancies), E""^^, plus 
the total energy of an atom at the kink site minus the 
total energy of the clean and unreconstructed surface, 

£,vac ^ ^vac ^ ^^Al-atom _ ^^^^^^ _ ^Al (111) ^^g) 

The vacancy formation energy has typically a positive 
value, as it usually costs energy to create a vacancy. In 
the second process, the vacancies are occupied by the 
adatoms. This allows us to define the binding energy per 
substitutionally adsorbed adatom, E^^^, 

Ijisub j-,Na/Al (lll)-sub j^ivac 

H = ^ad +-^/ ■ (1') 

For the on-surface adsorption the binding and adsorption 
energies are the same. 

In general, when a surface structure is considered 
which has a different number of substrate atoms to the 
clean unreconstructed surface the required or surplus 
substrate atoms are taken from or to kink sites at steps, 
using again the fact that in thermal equilibrium kink 
sites represent the thermodynamic reservoir for substrate 
atoms and establish that the chemical potential equals 
the cohesive energy. 

1, Atomic geometries of K on Al(lll) forO < Ok < 0.33 

This subsection describes various adsorbate geometries 
which are assumed by potassium atoms on the (111) sur- 
face of aluminum, depending on the adsorption tempera- 
ture and the adsorbate coverage. At first, non-activated 
adsorption is discussed, for which the calculations pre- 
dict that at low coverage the adatoms form a disordered 
homogeneous adlayer with the potassium atoms occupy- 
ing preferentially on-surface hollow sites. At 0k ~ 0.15 a 
structural phase transition occurs from this homogeneous 
adlayer into close-packed islands in which the adatoms 
occupy on-top positions. If thermal energy is available 
to overcome energy barriers, the adatoms in these is- 
lands undergo a site change into a substitutional geome- 
try. Surprisingly, the surface periodicity is unchanged in 
this major disruption of the surface, in which one third of 
the surface Al atoms are removed. The subsection is con- 
cluded with a brief comparison of the DFT-LDA results 
with recent experimental studies. 

The calculated adsorption energy for K on Al ( 1 1 1 ) for 
several conceivable adsorption sites and for various cov- 
erages are collected in Figs. ^ and where the latter 
also gives the binding energies for substitutional adsorp- 
tion, El^^^, and the vacancy formation energies, EJ^'^. 
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These results show that at low coverages the on-surface 
threefold hollow position is energetically preferred. The 
rapid decrease of the adsorption energy with increas- 
ing coverage reflects a strong adsorbate-adsorbate repul- 
sion. This is consistent with the strong adsorbate in- 
duced dipole moments which repel each other. The K 
atoms are expected therefore to be homogeneously dis- 
tributed over the surface at low coverages. For higher 
coverages the results indicate that a phase transition oc- 
curs to adsorbate islands with a condensed structure of 
(•\/3 X •\/3)i?30° periodicity. The local coverage in these 
islands is Qk = 1/3. Interestingly, the adsorption en- 
ergies for this periodicity are more favorable than those 
with a (2 X 2) periodicity and coverage 6k = 1/4. It 
can be therefore concluded that for coverages smaller 
than 0K ~ 0.15 the Gurney picture and the dipole- 
dipole repulsion (and depolarization) appear to describe 
the physics of the adlayer appropriately. However, when 
the coverage is increased, there is a crossover point, where 
a different mechanism becomes energetically preferable: 
Instead of keeping the largely ionic adsorbate-substrate 
bonds and the energetically unfavorable repulsive inter- 
action between adatoms, it becomes energetically ben- 
eficial to build up a close-packed potassium layer with 
a metallic bonding. Together with this, a reduction of 
the strength of the adsorbate-substrate bonding takes 
place. We had noted above that this is a very general 
phenomenon (roughly the bond strength scales with the 
square root of the coordination number). Thus, with the 
increase of coordination, the adsorbate-substrate bond 
becomes weaker and longer. The calculations show that 
in this condensed structure with (-s/S x \/3) periodicity, 
the on-top position is now energetically slightly more fa- 
vorable than the on-surface hollow position. Thus, for 
non-activated adsorption the theory predicts a structural 
phase transition from a homogenous, rather open adlayer 
with adatoms preferentially occupying the on-surface hol- 
low sites, into two-dimensional islands with a (\/3 x a/S) 
structure where the adatoms occupy on-top sites. This 
phase transition is indicated by the dashed line in Fig. |^a. 

Figure ^ shows that the adsorption energies of all the 
geometries considered with the (-s/S x \/3)i?30° period- 
icity are very similar, which indicates that at room tem- 
perature substitutional adsorption becomes energetically 
possible. In order to reach the substitutional geometry 
an energy barrier must be overcome to reconstruct the 
substrate surface. We will at first analyze the on-surface 
(%/3 X V3)i?30° adlayer and then we come back to the 
substitutional adsorption. 

For low temperature (i.e., non-activated) adsorption. 
Fig. ||a predicts that simultaneously with the island for- 
mation a site change occurs; the adatoms switch from on- 
surface hollow sites to on-top sites, which also involves a 
noticeable relaxation of the substrate H^,^ . In particu- 
lar, the Al atoms in the first substrate layer beneath the 
on-top K atoms are moved down and the Al atoms, be- 
tween the on-top K adatoms are shifted upwards, relative 
to the geometry of the unrelaxed substrate surface (see 



Fig. [7|a). It is important to note that the atomic radius 
of potassium is significantly larger than that of Al. This 
suggests, what is fully supported by the calculations, that 
K adatoms experience a rather small corrugation of the 
Al (111) surface. In fact, with the above-mentioned sub- 
strate relaxation, and because of the smaller K-Al bond 
length in the on-top geometry (due to the lower coordi- 
nation of the K adatom) potassium at on-top sites come 
even closer to the surface than for the on-surface hol- 
low sites. As a consequence, a better embedding of the 
adatoms in the substrate electron density is achieved. 

In the substitutional geometry, which can be reached 
when thermal energy enables energy barriers to be over- 
come, it is found that because of the larger size of the 
alkali metal atoms compared to the Al atoms, the K 
adatoms jut out of the surface (see Fig. |^). The effective 
radius of K in the substitutional site is larger than it is 
in the on-top site as can be seen clearly from comparing 
Figs. ^ and Qb. This trend reflects again the relation be- 
tween the bond strength and the local coordination: the 
higher the coordination, the weaker each bond is, and the 
longer each bond length is. The DFT-LDA results |2|] 
gives this increase in bond length as about 16%, and a 
LEED analysis Q gives 19%. With respect to the in- 
teratomic distances of the on-top geometry, the nearest- 
neighbor K-Al distances are 3.38 and 3.23 A, and for the 
substitutional site the values are 3.70 and 3.58 A, as ob- 
tained from the DFT-LDA calculations and the LEED 
analysis, respectively. A rough theoretical estimate of 
the temperature at which the transition from on-top to 
substitutional adsorption could occur with a measurable 
rate was obtained |32(| by assuming that the transition 
state of the reaction path is such that the K atoms are 
in substitutional sites and the ejected Al atoms are in 
isolated positions on the surface (not yet rebounded at 
steps). This transition state corresponds to an energy 
barrier of E'' = 0.8 eV. The Arrhenius equation, with 
a reasonable (but rough) estimate of the pre-exponential 
factor for the kick out event, then yields a transition tem- 
perature of 220 K. 

The nature of the adsorbate-adsorbate interaction for 
the substitutional site and the mechanism actuating the 
island formation with a substrate reconstruction with the 
(\/3 X \/3)_R30° structure can be learnt from Fig. 
It shows that the vacancy formation energy is particu- 
larly low for the x ^/3)R30° surface vacancy struc- 
ture compared to lower concentration vacancy structures. 
This behavior reflects an attractive vacancy-vacancy in- 
teraction, which is actuated by the fact that the group 
III element Al can easily assume covalent bonds. For 
the X \/3)-R30° surface vacancy structure, which is 
equivalent to that of a graphite layer, the coordination 
of the surface Al adatoms indeed gives rise to the forma- 
tion of directional bonds between the surface Al atoms. 
Figure ^ shows that the binding energy of K in a sur- 
face vacancy is relatively insensitive to the coverage and 
even exhibits a slight repulsive interaction between the 
K adatoms as indicated by the small decrease in binding 
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energy of the (\/3 x \/3)i?30° structure as compared to 
the lower coverage structures. Thus it is the relative sta- 
bility of the X v^)i?30° surface vacancy structure, 
which is the driving force for the island formation with a 
substitutional (VS x V3)i?30° structure 

Turning now to experimental studies we like to men- 
tion the low energy electron diffraction (LEED) study 
p2| , of which some results were already noted above. 
This investigation clearly proved that the (-y/S x \/3)R30° 
structure formed at low temperature contains K atoms in 
the on-top site, which, on warming to room temperature 
transforms irreversibly into a (V3 x \/3)R30° structure 
where the K atoms adsorb substitutionally. That is, even 
if the temperature is lowered the K atoms will remain in 
substitutional sites. 

The two different phase transitions predicted by the 
theory (on-surface hollow sites on-top site island 
formation, and on-surface — > substitutional geometry) 
are also consistent with high resolution core level spec- 
troscopy (HRCLS) measurements which show that 
island formation starts at Qk = 0.12 and 0.10 at low 
and room temperature, respectively. Similarly, a recent 
experimental study of the temperature dependence of 
optical second-harmonic generation |7^] found that at 
coverage Qk =0.06 no phase transition occurred when 
warming to room temperature whereas at 0k = 0.12 the 
transition took place. These experiments thus confirm 
the theoretical estimate, which was Qk ~ 0.15, quite 
well. The temperature at which the phase transition into 
the substitutional adsorbate islands takes place was mea- 
sured as 210 K which is in s urp rising agreement with 
the rough theoretical estimate noted above. 

2. Atomic geometries of Na on Al(lll) for < 0Na < 0.56 

For Na on Al (111) the same adsorbate geometries and 
coverages have been investigated by DFT-LDA calcula- 
tions as for K on Al(lll), with the addition of a struc- 
ture with a (4 X 4) periodicity corresponding to coverage 
©Na = 9/16 which represents a densely packed hexago- 
nal adlayer with nine Na adatoms per surface unit cell 
all occupying different (mostly low symmetry) on-surface 
sites. Furthermore, a variety of many different structures 
at ©Na = 0.5 was analyzed, because experiments had in- 
dicated a strongly intermixed surface geometry, without, 
however, giving a clue about the structure. The calcula- 
tions predict phase transitions which are very similar to 
those of K/Al (111), but some aspects are noticeably dif- 
ferent, due to the fact that Na is smaller than potassium 
although still larger than Al. The 0Na = 0.5 structure 
represents the most significant difference to the potas- 
sium adsorbate, but it is also unique with respect to 
any other adsorbate system: This structure may be de- 
scribed as an ordered four-layer surface alloy on an un- 
reconstructed Al(lll) substrate. 

We start with a discussion of the theoretical results 



for the non-activated adsorption; they are summarized in 
Fig. ||a and can be described as follows: For low coverages 
the on-surface hollow site is energetically most favorable 
and, in a manner similar to that for K on Al(lll), a 
strong repulsive adsorbate-adsorbate interaction is obvi- 
ous. In this case, as opposed to K on Al(lll) the on- 
top site is significantly less favorable than the on-surface 
hollow site at all coverages considered. This can be at- 
tributed to the smaller size of the Na atoms which expe- 
rience a more highly corrugated potential energy surface 
due to the substrate atomic structure than the larger 
K atoms. Figure ||a shows that the condensed struc- 
ture (the atomic geometry is displayed in Fig. 9a) for 
©Na = 9/16 is energetically more favorable than homo- 
geneous adlayers of Na for coverages larger than approx- 
imately ©Na = 0.1. Accordingly, these results imply 
that for very low coverages the adsorbates occupy on- 
surface hollow sites and are uniformly distributed over 
the surface (homogeneous adlayers) but for coverages 
greater than about ©Na = 0.1, island formation with 
the condensed structure and a (4 x 4) periodicity occurs 
as indicated by the dashed line. This is explained in 
much the same way as was done in the discussion of the 
condensation-into-islands phenomenon in Section IV A 1 
for K/A1(111). 

The theoretical results corresponding to activated ad- 
sorption are collected in Fig. Hb where for Na in the 
substitutional geometry the adsorption, binding, and va- 
cancy formation energies are displayed. It can be seen 
that the adsorption energy for the substitutional ge- 
ometry is the most favorable for all coverages investi- 
gated. In particular, substitutional adsorption with a 
(•\/3 X •\/3)i?30° periodicity has the most favorable ad- 
sorption energy. It is therefore expected that a conden- 
sation into (-s/S X V3)i?30° islands with the Na atoms in 
substitutional sites occurs, beginning at very low cover- 
ages, provided that thermal energy enables the transition 
to occur. The atomic geometry is depicted in Fig. 9b; it 
is equivalent to that of Fig. 7b with the only difference 
being the sizes of the adatoms. The mechanism which 
stabilizes the (a/S x V3)-R30° structure of the substitu- 
tional adsorption is identified as being the same as that 
for K on Al (111) [0, i.e., as being due to the particular 
stability of the (^/3x V3)i?30° surface vacancy structure. 

Upon further deposition of Na at room temperature 
onto the (VS x V3)i?30° substitutional adsorbate struc- 
ture, the formation of a (2 x 2) surface structure has 
been observed experimentally. This phase corresponds 
to a coverage of ©Na = 1/2 where there are two Na 
atoms in the surface unit cell. It has been the subject 
of a number of studies, and many different models have 
been proposed. Initially it was suggested that the struc- 
ture corresponds to three domains of (2 x 1) periodicity, 
involving Na atoms adsorbed on an unreconstructed sub- 
strate |73] . Subsequently it was argued that the structure 
consists of two layers of Na atoms, each with (2 x 2) peri- 
odicity, located on an unreconstructed substrate Efl] . Re- 
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cent studies using HRCLS |4^, normal incidence stand- 
ing x-ray wave-field (NISXW) and surface extended x-ray 
adsorption fine structure (SEXAFS) have surpris- 
ingly shown that a strong intermixing of the Na and Al 
atoms in the surface region occurs. In the NISXW study 
a model was proposed involving two reconstructed layers 
each of stoichiometry NaAl2 while from a recent scanning 
tunnel microscopy (STM) study [^4|j7|] it was suggested 
that the (2 x 2) phase consists of two Na layers where one 
Na adatom in the surface unit cell is located in a surface 
substitutional site and the second Na atom occupies an 
on-surface hollow site. In order to clarify the situation 
regarding the exact surface atomic structure of this con- 
troversial phase, DFT-LDA calculations were carried out 
for a large number of conceivable geometries . These 
calculations provided evidence for the uniqueness of the 
structure and permitted the rejection of all of the above- 
mentioned suggested models. The identified structure is 
significantly more complex and unexpected than any of 
the alkali metal adsorption systems studied so far. It is 
drawn in Fig. ^ and may be described as a composite 
double layer ordered surface alloy and has the lowest to- 
tal energy. In this structure an Al atom is incorporated 
into a double Na layer of which one of the Na atoms in 
the surface unit cell occupies a substitutional site and the 
second Na atom plus an Al atom are situated in fee- and 
hep-hollow sites, respectively. Alternatively, it could be 
described as a four layer Al-Na-Al-Na alloy on an unre- 
constructed Al(lll) surface. 

Intuitively one might expect that the Al atom in the 
composite double layer would adsorb in the fcc-hoUow 
site rather than in the hcp-hoUow site (continuing growth 
of the bulk), and that the Na atom would assume the hep- 
hollow site. It has been shown from the calculations that 
the reason this does not occur is due to a larger relaxation 
of the reconstructed Al layer (the Al layer containing 
three Al atoms per surface unit cell) and of the substrate, 
compared to the structure where the atoms occupy the 
opposite hollow sites. These relaxations correspond to a 
significant energy gain of 0.12 eV per surface unit cell. 

The adsorption energy per Na atom in this surface 
alloy structure is less than that of the {^/3 x \/3)i?30° 
substitutional phase and greater than that of the non- 
activated (4 X 4) surface geometry. The fact that it 
is greater than the (4 x 4) structure explains why the 
(4 X 4) structure is not found at room temperature. 
The result that the (2 x 2) surface alloy structure has 
a smaller adsorption energy than the substitutional ad- 
sorbate {VSx \/3)i?30° structure, implies that the (2 x 2) 
surface alloy will only form when it is not possible for the 
(73 x V3)mo° structure to grow any further. 

From consideration of the atomic structure of the 
(2 X 2) phase, it would seem that no mass transport is 
necessary in its formation. However, the lower cover- 
age (-\/3 X -\/3)_R30° substitutional structure involves dis- 
placement of 1/3 of a monolayer of Al atoms, which are 
assumed to diffuse across the surface to be re-adsorbed 
at steps. That is, the steps act as sinks for the displaced 



substrate atoms. The results indicate, therefore, that for- 
mation of the (2 X 2) structure from the (VS x \/3)i?30° 
structure involves the reverse process, that is, diffusion 
of 1/3 of a monolayer of Al atoms back from the steps 
which are used in the formation of the (2 x 2) structure. 
In this case the steps serve as sources of substrate atoms 

H- 

We note that the direct involvement of steps is one, but 
not the only possibility. Another mechanism could be 
that the kicked out Na atoms diffuse to regions between 
on-surface hollow-site Na atoms. This mechanism gives 
rise simultaneously to two domains of substitutional ad- 
sorption, one on the lower terrace from where the Al atom 
was ejected and one on the upper terrace to where the 
ejected Al atom finally resides. This two-domain mech- 
anism has been discussed for the Na adsorption on the 
Al(lOO) surface ||. 

It is appropriate to briefly mention some of the recent 
experimental studies. The results of the calculations for 
non-activated adsorption are compatible with a SEXAFS 
study of the (4 x 4) structure formed at low temperature 
which found that the Na atoms form a quasi-hexagonal 
close-packed adlayer where, within the unit cell, there is 
a slight clustering and rotation of the Na adatoms with 
resp ect to a strictly hexagonal close-packed arrangement 
|l76| . In the calculations these small displacements of the 
Na atoms within the (4 x 4) hexagonal adlayer were not 
considered. For activated adsorption the results of the 
DFT-LDA calculations are in full accordance with SEX- 
AFS |§, LEED H, and NISXW |^ studies of the 
(\/3 X \/3)R30° phase. Nearest neighbor Na-Al distances 
in the substitutional geometry of 3.13, 3.21, 3.31, and 
3.09 A were determined from the DFT-LDA calculations, 
LEED intensity, SEXAFS, and NISXW studies, respec- 
tively. The theoretical results are also in accordance with 
STM pictures taken at room temperature in the coverage 
range below 0Na = 1/2, which showed that islands with 
a (VS X •\/3)i?30° structure are present at low coverages 
and increase in size with coverage up to ONa = 1/3. The 
ordering of the total energies of the (^3 X V3)i?30° sub- 
stitutional structure and the (2 x 2) four- layer surface 
alloy structure is consistent with an STM study [ j74|j75[ 
which shows that initially islands with a (a/3 x \/3)R30° 
structure form rather than with a (2 x 2) structure. 

The complex structure of the four-layer surface alloy 
proposed in the DFT-LDA study has been confirmed by 
a detailed LEED intensity analysis SEXAFS data 

are consistent with the surface structure but the analysis 
does not allow a unique discrimination against alterna- 
tive models |52 . In addition, results of an x-ray photo- 
electron diffraction (XPD) study also agree with the 
predicted structure. For the first three interlayer spacings 
the values determined by DFT-LDA, LEED, SEXAFS, 
and XPD agree to within 0.13 A Q. The DFT-LDA 
and LEED results show, in particular, remarkable agree- 
ment. Both identify the presence of small (« 0.04 A ) 
lateral displacements of the Al atoms in layers four and 



10 



five and small vertical displacements of the Al atoms in 
layers five and six (counting the uppermost Na layer as 
the first, see Fig. |9p). The structure derived from the 
DFT-LDA total energy minimization is at variance with 
that suggested by the STM analysis of Brune et al. [ [74|j75| 
which is apparently due to the fact that STM is unable to 
detect all surface atoms with equal clarity and to identify 
the atomic species unambiguously. 

3. Atomic geometries of Na on Al (001 ) for < 0Na < 0.5 

The atomic density of the fee (001) surface is about 
15% lower than that of the (111) surface, and the (001) 
surface therefore has higher diffusion energy barriers to 
adatoms. As a consequence it is to be expected that the 
behavior of adsorbates on the (001) surface may exhibit 
some differences to on the (111) surface. We will show 
and explain in this subsection that at low temperatures 
Na adatoms always occupy on-surface hollow sites with 
a repulsive intra layer interaction and no indication of a 
phase transition (at low coverage) into close-packed is- 
lands. At higher temperature we find that Na adatoms 
occupy at low coverage also on-surface hollow sites, but 
at ONa ~ 0.15 a transition into islands with a substitu- 
tional geometry and a c(2 x 2) periodicity is predicted. In 
contrast to the other sections, we will discuss at greater 
length the results of recent experiments because the ex- 
perimental situation and published conclusions are more 
complex and contradictory than those for the other sys- 
tems. 

DFT-LDA calculations were carried out for many pos- 
sible surface geometries and coverages |3^,^ and some 
of these results are shown in Fig. [lC|a. 

We start with a discussion of the on-surface geometries 
which are relevant for non-activated adsorption. The en- 
ergies of different on-surface sites differ more widely than 
those for Na on Al (111) which can be attributed to the 
more corrugated potential energy surface of Al (001). For 
the c(2 X 2) structure (8Na = 0.5) the adsorption ener- 
gies for the on-top, bridge, and on-surface hollow sites 
are 1.33, 1.45, and 1.59 eV |36|. Therefore, only the on- 
surface hollow sites were considered further, the other 
sites were regarded as unlikely in the considered cover- 
age range. Figure [To| a shows that the on-surface hollow 
site is clearly preferred over the substitutional site at low 
coverages. Its adsorption energy rapidly decreases with 
increasing coverage, similar to K and Na on Al (111), in- 
dicating a strong repulsive interaction between the Na 
atoms. However, differently to the adsorption on the 
(111) surface, the theoretical results for Na on Al(OOl) 
give no indication that a condensation into close-packed 
islands takes place. This is quite certain for low cover- 
age, but for 8 > 0.3 we would not rule out that a phase 
transition into 8 = 0.5 condensed islands can occur. Fig- 
ure [lO| a suggests that the adlayer is compressed gradu- 
ally up to a coverage of 8Na = 0.5 for which Fig. |Tl|a 



displays the structure. The difference to the Al(lll) 
surface arises because the adsorbate-substrate bonding is 
stronger on Al (001) due to the better embedding of the 
Na atom when occupying the more open on-surface hol- 
low site. It therefore remains energetically advantageous 
to keep those strong bonds with the substrate rather than 
to build up metallic bonds in a condensed adlayer island. 

Figure [lO| a also contains the results relevant for acti- 
vated adsorption, namely the adsorption energies of sub- 
stitutional geometries. It can be seen that the adsorp- 
tion energy for Na in the surface substitutional site de- 
pends much more weakly on coverage than that for the 
on-surface Na, and in fact, the adsorption energy for the 
c(2 X 2) substitutional structure, which is displayed in 
Fig. [iTIb, is more favorable than that of the lower cov- 
erage substitutional structures. These results are inter- 
preted as follows: At low coverages the on-surface ad- 
sorption in a homogeneous adlayer is the stable structure 
for low as well as high temperature, but for higher cover- 
ages on-surface adsorption becomes metastable. For high 
temperature adsorption the adatoms then switch to sub- 
stitutional sites, forming islands with a c(2 x 2) structure 
(see Fig. |ll|b) which increase in size with increasing cov- 
erage. The phase transition from the on-surface hollow 
site to substitutional adsorption is indicated in Fig. p^ a 
by the dashed line. 

The shape of the adsorption energy versus coverage 
curve for Na in the substitutional site {E^^ in Fig. ^^) 
appears, at first glance, somewhat similar to the result 
for K and Na on Al (111) (see Figs, and^). However, 
from a plot of the adsorption, binding, and vacancy for- 
mation energies, which are also shown Fig. ^0|d, it can be 
seen that even though the adsorption energy versus cov- 
erage curve exhibits the same trend as for K and Na on 
Al (111), that is, in both cases the most favorable adsorp- 
tion energy is obtained for the highest coverage structure, 
the constituent components, namely, the binding and va- 
cancy formation energies comply to a different behavior. 
We see that on the (001) surface the binding energy is 
particularly favorable for the c(2 x 2) structure, reflecting 
an attractive adsorbate-adsorbate interaction which actu- 
ates island formation. It can be seen in addition that the 
surface vacancy formation energy shows little dependence 
on surface vacancy concentration. 

An additional surface structure was considered in the 
DFT-LDA calculations, corresponding to coverage 8 = 
0.5 which contains Na atoms in on-surface hollow and in 
substitutional sites js^]. The surface unit cell used was 
twice that of the c(2 x 2) cefl, i.e., (V2 x 2V2)i?45°. In 
this arrangement the Na atoms form a quasi-hexagonal 
structure, which could also be described as alternating 
strips of Na atoms in on-surface hollow and substitutional 
sites. It was found that the adsorption energy of this 
structure is between that of the on-surface hollow and 
substitutional sites. Accordingly, similar mixed hoUow- 
substitutional geometries may well be present in the tem- 
perature and coverage ranges at which the phase transi- 
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tion takes place, and they may also exist as metastable 
surface geometries. 

We now turn to a brief discussion of experimental stud- 
ies. Early LEED intensity analyses [||,|o) of c(2 x 2) 
Na/Al(001) concluded that the adatoms occupy the on- 
surface hollow site. In the preparation of the surface from 
which the intensity versus energy curves were measured 
the sample was annealed at 360 K. Recent experimen- 
tal results are, however, at variance with this conclusion. 
HRCLS |l| and SEXAFS || studies demonstrated con- 
vincingly that the atomic structure of the c(2 x 2) phases 
are different when formed at low temperature and at 
room temperature, that the structural phase transition is 
irreversible, and that the room temperature structure is 
due to a mixing of Na and Al atoms. The SEXAFS study 
determined the on-surface hollow site as the adsorption 
site at low temperature and for the structure formed at 
room temperature three possible models were proposed 
with one being favored which consists of Na atoms in 
on-top sites beneath a c(2 x 2) Al layer where the Al 
atoms are adsorbed in hollow sites with respect to the Na 
adlayer. The DFT-LDA analysis of these three surface 
alloy structures showed that they are energetically very 
unfavorable, and these structures have therefore been in- 
dubitably rejected [Q. Thus, these experimental results 
for the low temperature adsorption agree with the DFT- 
LDA results of Fig. |l^a. Also the experimental evidence 
for a phase transition towards an intermixed structure 
agrees with the theoretical prediction. However, the ge- 
ometry which was derived from the SEXAFS analysis 
is at clear variance with the DFT-LDA result shown in 
Fig. |ll|b. Interestingly, in an XPD study @0 of the 
room temperature c(2 x 2) phase of Na on Al (001) it was 
concluded that the surface may contain two domains, one 
with Na atoms in the substitutional site (45% occupancy) 
and the other with the Na atoms still in the on-surface 
hollow sites (55% occupancy). 

A recent dynamical LEED intensity analysis shows 
that the low temperature c(2 x 2) phase contains Na 
in on-surface hollow sites and in the room temperature 
c(2 X 2) phase the Na atoms adsorb substitutionally, in 
very good agreement with the DFT-LDA calculations. 
Importantly, the LEED study reveals that the transi- 
tion from on-surface hollow sites to substitutional sites 
proceeds in a gradual step-by-step way, beginning at a 
temperature of approximately 180 K, and that on com- 
pletion of the phase transition (determined to occur at 
about 260 K) the contribution to the LEED intensities 
from substitutional adsorption is not less than 90%. 

As discussed in Ref. one particular atomic config- 
uration consistent with the behavior of the temperature 
dependent satellite features, which are seen during the 
phase transition, is one which corresponds to the cre- 
ation of domains arranged in narrow rows or "strips" 
where the Na atoms are in substitutional sites. As the 
temperature is increased and the phase transition pro- 
ceeds, the rows or strips of Na atoms in substitutional 
sites become progressively wider until at the end of the 



transition all adatoms are in substitutional sites, recov- 
ering again the sharp c(2 x 2) LEED pattern |Q. Three 
possible mechanisms are envisaged by which the substi- 
tutional geometry could be realized: i) The Na atoms 
cause the ejection of substrate Al atoms which diffuse 
across the surface to steps - the steps acting as sinks. 

ii) Steps on the surface act as sources of Al atoms which 
diffuse to the regions between the on-surface hollow-site 
Na atoms, forming exactly the substitutional geometry. 

iii) As in i) where the Na atoms eject substrate Al atoms 
which, however, do not move to steps but diffuse as in ii) 
to the regions between on-surface hollow-site Na atoms. 
The latter mechanism gives rise simultaneously to two 
regions of substitutional adsorption, one on the "lower" 
terrace from where the Al atom was ejected and one on 
an "upper" terrace to where the ejected Al atom finally 
resides. "Kick-out" of the Al atoms by the Na atoms 
could proceed in a "chain-reaction-like" process which 
would give rise to the formation of rows or strips of Na 
atoms in substitutional sites, either on one terrace only if 
process i) occurs or on two terraces if process iii) occurs. 
The starting point for the "kick-out" process could con- 
ceivably occur at imperfections or steps on the surface. 

The interatomic distances obtained by the various ex- 
periments are in good agreement with those obtained 
from the DFT-LDA calculations. For the on-surface 
hollow site, the obtained values read 3.08, 3.21, 2.93, 
3.27 A as obtained by the DFT-LDA, SEXAFS, XPD, 
and LEED studies. For the room temperature phase the 
Na-Al bond lengths determined by the DFT-LDA, XPD, 
and LEED investigations are 3.02, 3.11, and 3.07 A re- 
spectively. For the substitutional site, DFT-LDA and 
LEED identify a small (« 0.05 A ) vertical displacement 
of the Al atoms in the second Al layer. Both studies 
also find that the effective radius of Na in the substitu- 
tional site is shorter than in the on-surface hollow site: 
approximately 4% from DFT-LDA and 11% from LEED. 

In a manner similar to that described above for K on 
Al, the temperature for the on-surface hollow-to-surface 
substitutional transition has been estimated from DFT- 
LDA calculations. The value obtained was 180 K and 
calculated by assuming that the transition state of the 
reaction path is such that the Na atoms are in substitu- 
tional sites and the ejected Al atoms are in the energet- 
ically least favorable position of the exchange diffusion 
mechanism (two Al atoms above a surface vacancy, see 
Refs. ^,Q), not yet rebounded at kink sites at steps. 
This rough theoretical estimate of the transition tem- 
perature is in good agreement with the value of 160 K 
determined from the HRCLS stud y |8l[ and with 180 K 
as determined in the LEED study 



B. Adsorbate bonding and electronic properties 

The widely differing adsorbate geometries and their 
strong dependencies on coverage and temperature, as de- 
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scribed in the previous section, go together with signif- 
icant changes of the surface chemical bond and surface 
electronic structure. In this respect, in this section, we 
discuss the adsorbate-substrate and adsorbate-adsorbate 
interaction of the different surface structures, the associ- 
ated electron densities at the surface, and the resulting 
work function changes. 

As mentioned in the introduction, the characteristic 
change in the work function with coverage, A<I>(0), is 
typically explained in terms of the Gurney picture and 
the "usual form" (see Fig. ^ can be described reasonably 
well within a model in which a high Ts jellium is adsorbed 
on = 2 jellium. Alternatively, but with essentially 
the same conclusions, Muscat and Newns 1^,^ exam- 
ined the coverage dependence of the adsorbate-substrate 
and adsorbate-adsorbate interaction using an empirical 
Anderson- type Hamiltonian js^,^^, which takes the al- 
kali metal valence states and the eigenstates of the semi- 
infinite substrate properly into account. They studied 
the work function change for alkali metal atoms on tran- 
sition metal surfaces . 

To obtain quantitative agreement between calculated 
and experimental results it is necessary to take into ac- 
count in the calculations that the atomic geometry (of 
adsorbate and substrate) and the surface electronic struc- 
ture are closely linked. Thus, a realistic description of 
the surface electronic properties necessarily requires a 
correct treatment of the various adsorbate arrangements 
discussed in Section IV A above. To highlight this point 
we investigate the effect on the calculated adsorbate in- 
duced change of the work function by assuming hypothet- 
ical atomic arrangements of Na on Al (111) where the Na 
adatoms are homogeneously distributed over the surface 
for the specified coverage. The self-consistently calcu- 
lated surface dipole moments and work function change 
corresponding to these structures are shown in Figs, [l^ a 
and 123, respectively. In agreement with the traditional 
picture of alkali metal adsorption, there is a significant 
decrease of the surface dipole moment with increasing 
coverage. For comparison, the results for homogeneous 
adlayers of Na in the substitutional site are also shown. 
Here we see that the depolarization effect, which is one 
indicator of the strength of the adatom-adatom interac- 
tion, is much less dramatic. This is due to the fact that in 
the substitutional site the Na atoms sit lower in the sur- 
face and the repulsive adsorbate-adsorbate interaction is 
screened better. As a consequence of the lower (closer to 
the substrate) position of the substitutional adsorbate we 
also see that the dipole moment (at low coverage) is much 
smaller - largely due to the smaller dipole length. As a 
consequence of this better screening of the adsorbate- 
adsorbate interaction for substitutional adsorption, i.e., 
the (near) absence of depolarization, the A$(8) curve 
for the homogeneous substitutionally adsorbed adlayer 
does not exhibit a minimum. The obtained work function 
change in Fig. [T^ agrees reasonably well, as expected, 
with that of Lang's jellium on jellium calculation. How- 
ever, as we will discuss below, the curves for Fig. and 



are in clear disagreement with the measured results. 



1. K on Al(lll) 

The calculated change in work function of potassium 
on Al (111) corresponding to the realistic atomic arrange- 
ments, which were discussed in Section IV A 1, are dis- 
played in Figs, and ^3|b for non-activated and acti- 
vated adsorption, respectively. For the non-activated ad- 
sorption the adatoms stay on the surface, but at 8 « 0.15 
they change from a homogeneous adlayer to close-packed 
islands (the structure was shown in Fig. 0a), which for 
increasing coverage increase in size. This phase transi- 
tion is the reason for the minimum in the dashed curve 
of Fig. p^. The comparison with the low temperature 
experimental result shows reasonable agreement. In par- 
ticular we emphasize that the origin of the minimum in 
A<I>(8) is not due to the adsorbate-adsorbate depolar- 
ization, but to the phase transition into close-packed is- 
lands. It is interesting to inspect the adsorbate- induced 
change of the electron density (see Fig. |lja). Within 
this structure, the adsorbate-adsorbate distance is only 
5% larger than that of bulk K, and Fig. [l4^ indeed re- 
veals that there is considerable charge density between 
the K adatoms. We therefore tend to classify the bond- 
ing within the adlayer as metallic. Consistent with this 
we note that the surface band structure of this adlayer is 
similar to that of a free-standing (\/3 x \/3)^30°-K layer 

The calculated work function change for activated ad- 
sorption is shown in Fig. ^3|b. Although this curve looks 
similar to that of the low temperature (non-activated) 
adsorption, this similarity is in fact highly incidental. 
The physics behind the curve of Fig. is that of a 
structural phase transition at 8 « 0.15 into islands with 
substitutional adatoms (compare Fig. ^p). Incidentally, 
this phase transition occurs at about the same coverage 
as the phase transition of the non-activated adsorption 
which implies that the minima of the A<I>(8) curves of 
Figs. [T^ a and 133 are at about the same positions. The 
charge density between the adsorbates (see Fig. p^), 
which is essentially provided by the substrate atoms, is 
significant and screens the adsorbate-adsorbate interac- 
tion quite efficiently. Therefore the substitutionally ad- 
sorbed potassium can attain its ionic character even in 
the close-packed (-s/S x -\/3)i?30° structure in clear con- 
trast to the on-surface adsorbate structures. The DOS 
for the substitutional adsorbate (not shown) is also sig- 
nificantly different to that of the on-surface adsorbate 
structure, and it is mainly due to the reconstruction of 
the Al substrate, i.e., to the removal of one third of all 
surface atoms; the K atoms which occupy the surface 
vacancies modify the vacancy electronic structure only 
little Wh- 
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2. Na on Al(lll) 

The behavior of Na on Al(lll) is quaUtatively simi- 
lar to that of K on Al(lll). However, because of the 
smaller size of the Na atom, the condensed structure of 
the non-activated adsorption is more dense for Na than 
for K. The (4 x 4) structure of the Na islands (local cov- 
erage = 9/16), which is assumed for low temperature 
adsorption and coverages above w 0.1 was displayed in 
Fig. The calculated work function change is shown 
in Fig. |l^ for different temperatures. The minimum at 
O « 0.1 marks the phase transition from the homoge- 
neous adlayer to the condensed islands which grow in 
size with increasing coverage. Again we emphasize that 
the shape of the A$(8) curve and the physics behind this 
shape is significantly different to that of the "traditional" 
picture of Figs. || and ^2|b. In the condensed structure the 
nearest-neighbor Na-Na distance is very similar to that 
in bulk Na and the intra layer bonding can be described 
as metallic. 

The calculated work function change with coverage for 
the activated adsorption follows a linear relation between 
A$ and (see Fig. [l5|b). This results, because at high 
temperature island formation starts at very low cover- 
age and thus the work function change is mainly deter- 
mined by the surface dipole moment of the condensed 
(-s/S X -\/3)i?30°-Na substitutional structure times the 
number of Na atoms in these islands. The work function 
decreases up to the point where the whole surface is cov- 
ered with (Vs X V3)-R30° substitutional Na atoms. For 
coverages higher than 8 = 1/3 (compare Section IV A 2) 
a major restructuring of the surface takes place. The 
analysis of the work function shows that A$(8) should 
exhibit a minimum at 8 = 1/3 (A<I> w —1.6 eV) and 
then increase up to A(f>(8) = -1.16 eV at 8 = 1/2 for 
the four-layer surface alloy. For even higher coverages 
the curve should reach the value of the work function of 
a Na crystal surface. 

Similar to that of K in the substitutional site, the bond- 
ing of Na in the (v^ x V3)^30° substitutional structure 
can be described as partially ionic. The density difference 
plot is shown in Fig. ^ 

On further deposition of one sixth of a monolayer of Na 
onto the substitutional (\/3 x V3)-R30°-Na surface, the 
(2 X 2) surface alloy structure forms. The adsorbate in- 
duced electron density is shown in Fig. |l^. Two regions of 
charge density maxima occur. Both are weighted towards 
the uppermost Al atom. Inspection of the total charge 
density reveals that the Al atom in the hcp-hoUow site 
forms covalent-like bonds with the three Al atoms under 
it. This is possible because the Al atoms in the recon- 
structed layer underneath are rather lowly coordinated 
due to the displacement of one quarter of a monolayer of 
Al atoms. The Al-Al bondlength is shorter than in bulk 
Al, compare 2.66 to 2.81 A (theoretical bulk value) re- 
flecting the stronger and lower coordinated bonding. The 
occupied adsorbate induced DOS for the (2 x 2) structure 



is largely determined by the Al structure, which corre- 
sponds to the Al(lll) surface with 8 — 1/4 vacancies 
plus 8 = 1/4 Al atoms in the on-surface hcp-hollow sites 
(compare Fig. j^c) . In particular we noted a peak at ap- 
proximately 2 eV below the Fermi level (see Ref. [Q). 

3. Na on Al(OOl) 

Considering now Na on Al (001), we see that from in- 
spection of the total charge density for the c(2 x 2) struc- 
ture with Na adsorbed substitutionally, that directional 
bonds are formed between Al atoms of the uppermost 
layer (from which 8 = 0.5 Al atoms have been removed) 
and their four nearest-neighbor atoms in the layer un- 
derneath. The density difference plot in Fig. ^ clearly 
displays a maximum of charge density located above the 
uppermost Al atoms. From inspection of the adsorbate 
induced DOS, an occupied and an unoccupied feature 
at about 2 eV below and 1 eV above the Fermi level, 
respectively are identified. While the peak below the 
Fermi level is largely due to the Al structure (substrate 
with vacancies), the peak above the Fermi level arises 
from the Na atoms. An analysis of the wave functions of 
these states indicated that (at least at F) the occupied 
state is localized above the uppermost Al atoms and the 
unoccupied state is localized above the Na atoms. 

V. CONCLUSION 

In the present paper we discussed results of recent den- 
sity functional theory calculations for alkali metal adsor- 
bates on close-packed metal surfaces. The analyses sup- 
port the traditional charge transfer picture of alkali metal 
adsorption as being an appropriate description in the 
limit of low coverage and non-activated adsorption, i.e., 
when isolated alkali metal adatoms occupy on-surface po- 
sitions. It is emphasized, however, that low coverage 
means very low because adsorbate-adsorbate interactions 
can become very noticeable already at surprisingly low 
coverages (e.g. 8 < 0.1) and the traditional picture, 
thus, should not be generally applied in the higher cover- 
age regime. In particular we stressed the interplay of the 
electronic properties of the surface, which is for example 
reflected in the work function, and the adsorbate site. 

The transition region from the regime at very low cov- 
erages, governed by electrostatic repulsion, to higher cov- 
erage regimes where adlayer phases develop which pos- 
sess a cohesive interaction between the adatoms (or va- 
cancies) is particularly interesting. In this respect, it 
has been shown that with increasing coverage, Na and 
K adsorbates on close-packed Al surfaces exhibit an un- 
expected variety of physical and chemical phenomena. 
For non-activated adsorption these phenomena include 
site changes and on-surface condensation phase transi- 
tions and for activated adsorption, condensation phase 



14 



transitions involving substitutional adsorption, as well 
as surface alloy formation. Such unusual and varied be- 
havior as described in the present work with respect to 
the knowledge of several years ago, emphasizes a new up- 
dated view of alkali metal adsorption on metal surfaces. 

The novel adsorption geometries were observed so far 
for the two close packed aluminum surfaces: Al (111) and 
Al(lOO). Thus, the metal, which so far was assumed 
to be the most jellium-like system, reveals in fact quite 
clearly the influence of its atomic structure and tendency 
towards geometries which are stabilized by the formation 
of partially covalent bonds. At this point it is not clear, 
whether similar behavior ought to be expected for other 
systems. Calculations of the formation energy of surface 
vacancies indicate that silver and maybe also copper 
are further candidates for substrates where substitutional 
adsorption might occur. In fact, substitutional adsorp- 
tion is a much more general phenomenon than hitherto 
expected. Atoms which are immiscible with the bulk 
may well intermix at the surface, because there the size 
mismatch can be easy accounted for by relaxations. 
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FIG. 1. Schematic representation of the emptying and 
broadening of the alkali metal atomic s level as it approaches a 
metal surface. $s and A"!" represent the substrate work func- 
tion and the change in the substrate work function, respec- 
tively, and Ef and la denote the Fermi level of the substrate 
and the first ionization energy of the free adatom. The figure 
is drawn with values which resemble the adsorption of Na on 
Al(lll) at low temperature and low coverage (O ~ 0.1). 



FIG. 2. Change in work function with increasing coverage 
obtained by Lang using a jellium on jellium model for 
parameters corresponding to Na on Al(lll). 
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FIG. 3. (a) Change in the density of states due to adsorp- 
tion of sodium (dashed), silicon (dot-dashed), and chlorine 
(dotted) on Al(lll) (from Ref. ||l[). The density of states 
of the aluminum substrate is shown as a solid line. Energies 
are given with respect to the Fermi level which is indicated 
by a vertical full line. The vertical dashed line indicates the 
bottom of the valence band, (b) Analogous quantities to (a) 
for a jellium substrate with an electron density which is ap- 
proximately equal to the average value of Al: rs = 2 bohr. Li 
is considered instead of Na (after Fig. 2 of Ref. j^). 

FIG. 4. Total valence electron density in a plane perpen- 
dicular to the surface in the [121] direction for the chemisorp- 
tion of sodium, silicon, and chlorine (from left to right) on 
Al(lll). The contour spacing is 29 x {l.^f x 10"^bohr"^ 
with —5 < k < +5. The green-yellow boarder line (fe = 0, 
n(r) = 29 X 10~^bohr~"^ is the average density of the alu- 
minum, i.e., Ts = 2.02 bohr. 



FIG. 5. Density difference, n(r) - n^ir) - n''*°'" (r) , for 
Na, Si, and CI on Cu(lll) (upper panel) and on Al(lll) 
(lower panel) in a plane perpendicular to the surface in the 
[121] direction (from Ref. ||). n(r), n°(r), and n''*°'"(r) 
are the electron densities of the adsorbate system, the clean 
surface, and the free, neutral adatom. The hatched areas 
correspond to regions of charge increase and the unhatched 
areas correspond to charge depletion. The units are 10""^ 
bohr~^ and the contour lines correspond to the values —5.0, 
-2.5, -0.5, 0.0, 0.5, 2.5, and 5.0 . 



FIG. 9. Perspective view of (a) the condensed hexagonal 
(4 X 4)-Na/Al(lll) surface structure, (b) (^ x V3)R30°- 
Na/Al (111) with Na in the surface substitutional site, and (c) 
the (2 X 2)-Na/Al (111) surface alloy structure. The purple, 
blue, and green circles represent, respectively, Na atoms and 
Al atoms in the first and second layer. In (c) the Al atom in 
the hcp-hoUow site of the surface alloy is shown in yellow and 
the inequivalent Na atoms are depicted by different shades of 
purple. 

FIG. 10. (a) Adsorption energy versus coverage for Na 
on Al (001) in the on-surface hollow site and in the surface 
substitutional site, (b) Adsorption, E^, and binding en- 
ergies, El^^, for substitutional adsorption, and the vacancy 
formation energy, versus coverage (from Ref. p7[). 

FIG. 11. Perspective view of the atomic structure of 
c(2 X 2)-Na on Al(OOl). The purple, blue, and green cir- 
cles represent, respectively, Na atoms and Al atoms in the 
first and second layer, (a) On-surface hollow site geometry 
and (b) substitutional geometry, where every second Al atom 
in the top unreconstructed Al layer has been kicked out. 

FIG. 12. (a) Surface dipole moment for homogeneous ad- 
layers of Na on Al(lll) in the fcc-hoUow site and in the sur- 
face substitutional site versus coverage, (b) The correspond- 
ing change in the work function where the same symbols as 
in (a) are used (after Ref. ]7C[]). 



FIG. 6. (a) Adsorption energy versus coverage for K in 
on-top and fcc-hoUow sites on Al(lll). The dashed line 
marks the phase transition from the homogeneous adlayer 
into adatom islands with a {y/S x \^)R30° periodicity (from 
Ref. ||^). (b) Adsorption, El'^^, and binding energies, El""^, 
for substitutional adsorption, and vacancy formation energy, 
-B™, versus coverage (from Ref. [Q ). 

FIG. 7. Perspective view of the atomic structure of 
(^x V3)i?30°-K on Al (111), (a) K atoms (red) are in on-top 
sites and the substrate atoms under the potassium atoms (yel- 
low) are displaced towards the bulk and the substrate atoms 
between the potassium atoms (blue) are displaced upwards 
with respect to the ideal positions, (b) Substitutional geom- 
etry where the Al atoms which were beneath the K atoms in 

(a) have been kicked out. 

FIG. 8. Adsorption energy versus coverage for Na in 
on-top and fcc-hoUow sites on Al(lll). The dashed line 
marks the phase transition from the homogeneous adlayer into 
adatom islands with a condensed structure (from Ref. [|3l|). 

(b) Adsorption, -Bid'', and binding energies, Er^, for sub- 
stitutional adsorption, and vacancy formation energy, EJ^, 
versus coverage (from Ref. ) . 



FIG. 13. Calculated (dashed line) and measured (full line 
with open circles) work function change, A$, with coverage 
for K on Al (111), (a) For "low temperature" and (b) for "high 
temperature" (from Ref. j3l]). The experimental results are 
from Ref. [||. 

FIG. 14. Change of the electron density for 

{VS X VS)RSO°-K on Al(lll) for K in (a) the on-top site 
and (b) the substitutional site. For the on-top site the ref- 
erence system is the clean surface and for the substitutional 
site it is the (-\/3 x \/^)R30° surface vacancy structure. The 
contours are in the plane normal to the surface and parallel 
to the [121] direction. The substrate atoms are represented 
by the small dots and the K atoms by the large dots. The 
units are 10"^ bohr"^ (from Ref. 

FIG. 15. Calculated (dashed line) and measured (full line 
with open circles) work function change, A$, with coverage 
for Na on Al(lll). (a) For "low temperature" and (b) for 
"high temperature" (from Ref. |3^). The experimental re- 
sults are from Ref. Ma. 



17 



FIG. 16. Change of the electron density for 

(VS X y3)7?30°-Na on Al(lll) for Na in the substitutional 
site. The reference system is the (y/S x y/S)R30° surface va- 
cancy structure. The contours are in the plane normal to 
the surface and parallel to the [121] direction. The substrate 
atoms are represented by the small dots and the Na atoms by 
the large dots. The units are 10"^ bohr~^ (from Ref. [|l|). 

FIG. 17. Change of the electron density for the 
(2 X 2)-Na/Al (111) surface alloy structure The reference sys- 
tem is the Al (111) surface plus the reconstructed Al layer and 
the Al atom in the hcp-hollow site on the reconstructed layer. 
The contours are in the plane normal to the surface and par- 
allel to the |TT2] direction. The Al atoms are represented by 
the small dots and the Na atoms by the large dots. The units 
are 10"^ bohr"^ (from Refs. Q). 

FIG. 18. Change 

of the electron density for c(2 x 2)-Na/Al (001) with the Na 
atoms in the substitutional site. The reference system is the 
c(2 X 2) surface vacancy structure. The contours are in the 
plane normal to the surface and parallel to the [110] direction. 
The Al atoms are represented by the small dots and the Na 
atoms by the large dots. The units are 10"'^ bohr^'' (from 
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